Topological Renyi Entropy After a Quantum Quench 
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We present an analytical study on the resilience of topological order after a quantum quench. The system is 
initially prepared in the ground state of the Toric Code Model, and then quenched by switching on an external 
magnetic field. During the subsequent time evolution, the variation in topological order is detected via the 
topological Renyi entropy of order 2. We consider two different quenches: the first one has an exact solution, 
while the second one requires perturbation theory. In both cases, we find that the long-term time average of the 
topological Renyi entropy in the thermodynamic limit is the same as its initial value. Based on our results, we 
argue that topological order is resilient against a wide range of quenches. 



Introduction. — Topologically ordered phases in quantum 
many-body systems are of extreme importance in condensed 
matter physics HI and in quantum information HH^t] . They are 
novel phases of matter that defy the Landau paradigm of spon- 
taneous symmetry breaking and possess a robust ground state 
degeneracy that makes them good candidates for quantum 
memories. Moreover, they feature anyonic excitations whose 
interactions are of a topological nature and that are therefore 
less subject to decoherence J2, 0]. Since these phases can- 
not be described by local order parameters, their characteriza- 
tion must be done using global properties of the ground state 
wavefunction. Recently, it has been argued that topological 
order in the ground state is revealed by an entanglement pat- 
tern called topological entropy: a universal correction to the 
boundary law for the entanglement entropy It has also 

been shown that this correction is robust against perturbations 
therefore it serves as an effective non-local order parameter 
for topologically ordered phases JHIil]. 

Unfortunately, the entanglement entropy is extremely hard 
to compute and its measurement requires a complete state to- 
mography of a possibly extended system Holl . On the other 
hand, there is now an increasing understanding of the fact that 
the Renyi entropy of order 2 contains an appreciable amount 
of information about the universal properties of a quantum 
many -body system Hill . In particular, the aforementioned 
topological pattern of entanglement appears in the Renyi en- 
tropy as well Jfj , 12 ] . Moreover, this quantity is significantly 
easier to compute and can in principle be measured directly 
1 1311 . Indeed, there are already proposals for the experimental 
realizations of such measurements 11411 . 

One of the most interesting questions regarding topological 
order is its behavior away from equilibrium. In general, the 
non-equilibrium properties of quantum many-body systems 
are among the most fruitful subjects in both theoretical 11511 
and experimental lfl6ll condensed matter physics, and they can 
be conveniently studied in the setting of the quantum quench: 
a sudden change in the parameters of the system Hamiltonian 
fl7ll . In particular, the quantum quench of a topologically or- 
dered system was numerically studied in Ref. Ilal . where they 
found that topological order is resilient against certain types of 
quenches. The main disadvantage of their method is that it is 
only applicable to small system sizes. 



In this Letter, we analytically study the behavior of a topo- 
logically ordered system after a quantum quench by deter- 
mining the time evolution of the topological Renyi entropy 
of order 2. In particular, we prepare the system in the ground 
state of the Toric Code Model (TCM) |:2j[] and quench it by 
switching on an external magnetic field. Concentrating on 
two different versions of the quench, we establish both an ex- 
act treatment and a perturbation theory. Due to the nature of 
these methods, our results are not confined to the small system 
sizes accessible by exact diagonalization H, 18]. 

General formalism. — We consider the TCM with an exter- 
nal magnetic field in the +z direction. In this model, there are 
2N 2 spins on the edges of an TV x N square lattice with peri- 
odic boundary conditions J2l- The spins on the horizontal (h) 
and the vertical (v) edges experience different magnetic fields 
A and kA, therefore the Hamiltonian takes the form 11211 



where the star operators A s = Yiies &f an< ^ tne plaquette op- 
erators Bp = Yii£ P &i belong to stars (s) and plaquettes (p) 
on the lattice containing four spins each (see Fig. [TJ. 



FIG. 1 : (Color online) Illustration of the square lattice with the phys- 
ical spins located at the horizontal (black circles) and the vertical 
(white circles) edges. Examples of a star (red cross labeled s) and a 
plaquette (blue square labeled p) are included. 

In the case of A = 0, the system is exactly solvable be- 
cause [A S ,B P ] = [A s ,H(0)} = [B p ,H{0)\ = 0. The sectors 
with different sets of expectation values {A s = (A s ) = ±1} 
and {Bp = (B p ) = ±1} can be treated independently. In 
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the ground-state sector with A s = +1 (Vs) and B p = +1 
(Vp), there are four linearly independent degenerate ground 
states that are distinguished by the topological quantum num- 
bers Z\ % = ±1. These quantum numbers are expectation 
values for products of erf operators along horizontal and ver- 
tical strings going round the lattice. The ground state with 
Z x = Z 2 = +1 becomes |0) = NJ[ S {1 + A s ))\ fr), where 
JV is a normalization constant and | ft ) is the completely po- 
larized state with erf = (erf) = +1 (Vi). 

In the case of A > 0, the system is not exactly solvable in 
general because [A s ,<rf] ^ 0. On the other hand, it is true 
that [B p ,Z xa ] = [B P ,H(X)] = [Z h2 ,H(X)] = 0, therefore 
the sectors with different values of B p and Z\p can be treated 
independently. In the following, we only consider the lowest- 
energy sector with B p = +1 (Vp) and Z\ = Z 2 = +1. The 
states within this sector are distinguished by the expectation 
values {A s = ±1}, and we introduce a corresponding repre- 
sentation with quasi-spins A s located at the stars |19J. In this 
quasi-spin representation, the operator A z s = A s measures 
and the operator A x s flips the quasi-spin A s . Up to an irrele- 
vant additive constant, the Hamiltonian in Eq. (Q~|i becomes 

H{\) = -Y,K-\ Y, Mai-kx a * A*" 



Fig. |2]as S% = -S ( a ] + S [ a 2) + S { a 3) - Sg>. When defining 
S^ in the standard way JsJ] , we take Sa = S AB between 
the two thick subsystems A and B. However, it was argued 
in Ref. Il2ll that the presence of the 1 2 lattice gauge struc- 
ture due to the constraint B p = +1 (Vp) makes it possible 
to substitute the subsystem A with its boundary dA. In the 
following, we use the corresponding modified definition for 
S% in which we take S a m) = A ' M = S® A between the 
thin boundary subsystem dA and the rest of the system. The 
topological Renyi entropy is non-zero if and only if the given 
state has topological order. For example, its value is S^ = 2 
for the TCM ground state |0) and S J = for the completely 
polarized state | f|- ) . Note that the dimensions D and d of the 
subsystems need to be macroscopic such that D > d ^> 1. 



(s,s')£h 



(2) 

where (s, s') means that the summation is over edges between 
neighboring stars s and s'. The TCM with external mag- 
netic field is therefore equivalent to a 2D transverse field Ising 
model (TFIM) in which the coupling strengths on the horizon- 
tal and the vertical edges are not equal in general. 

When studying the quantum quench, we are interested in 
the time evolution of the quantum state after a sudden 

change in the Hamiltonian. At time t = 0, the system is set up 
in the ground state of the initial Hamiltonian H(0) such that 
|*(0)) = |0). At time t > 0, the system is evolved with the 
quench Hamiltonian H (A) and the state takes the general form 
= exp[— itH (A)]|0). To extract any valuable informa- 
tion from this expression, we need to obtain a full solution 
of the Hamiltonian H(X). In the following, we consider two 
important limits. When k = 1, the horizontal and the verti- 
cal coupling strengths are equal, and the equivalent TFIM is 
the standard 2D TFIM. When k — >■ 0+, the vertical coupling 
strength vanishes, and the equivalent TFIM factorizes into N 
independent ID TFIM copies along the horizontal chains of 
the lattice I20I1 . The k = case has an exact solution available 
for all values of A, while the k = 1 case requires perturbation 
theory around the exactly solvable point at A = 0. 

Topological Renyi entropy. — The Renyi entropy of order a 
is a generalization of the usual (von Neumann) entanglement 
entropy that characterizes the quantum entanglement between 
two complementary subsystems A and B = A. It is defined 
by S AB = log 2 Tr[p2J/ (1 — a), where pa is the reduced den- 
sity operator for A. Note that the usual entanglement entropy 
is recovered in the special case of a = 1. 

The topological Renyi entropy is extracted from indepen- 
dent Renyi entropies Sa calculated in the four cases (m) of 
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FIG. 2: Illustration of the subsystems A and B in the four cases (m) 
that are used to calculate the topological Renyi entropy. 

The modified definition for the topological Renyi entropy 
is an immense simplification to our calculations because the 
reduced density matrix pa is then diagonal in the basis of the 
physical spins erf. It was shown in Ref. Ill 211 that if the bound- 
ary subsystem dA consists of n closed loops with a combined 
length L, the Renyi entropy of order 2 for a generic state | 
satisfying the gauge constraint B p = +1 (Vp) is 



S° A = 



i+Y.(*\a x s1 AiW 



(L-n)- log 2 x , ^ N 

L S1,S 2 

^\A%k^Aik^f 



E 

Sl ,52 ,83 ,S4 



(3) 



where the sum contains all 2 L ~ n possible products with an 
even number of quasi-spin operators A x s chosen from each 
closed loop of the subsystem dA. 

Exact treatment. — In the case of n — 0, there are no in- 
teractions between the horizontal chains of the lattice, and the 
2D system factorizes into N independent ID systems in terms 
of the quasi-spins A s . Indeed, the 2D Hamiltonian in Eq. (0 
becomes the direct sum of N identical ID Hamiltonians. If 
we consider any horizontal chain and label its N stars with 
1 < I < N, the corresponding ID Hamiltonian is 



h(X) 



N 

E 

1=1 



A? + XATA 



1 -"-i-i 



(4) 



where the periodic boundary conditions are taken into account 
by A = A N . On the other hand, the ID TFIM Hamilto- 
nian in Eq. can be solved exactly by following a standard 
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procedure 12111 , In the first step, we map the quasi-spins Ai 
to fermions via the Jordan— Wigner transformation given by 
1 - 2c|q and Af = (A+)t = c|e"^=i °h, where 
and q are standard fermionic creation and annihilation op- 
erators. The translational symmetry is then exploited by the 
Fourier transform q = N^ 1 / 2 J2k elkl(: k, where the sum is 
over the momenta k m = n(2m — l)/N. We finally intro- 
duce new fermionic operators via the Bogoliubov transforma- 
tion Cfc = cos 0k Jk + i sin k ^_ k , and the Hamiltonian in Eq. 

takes the free-fermion form h(X) = ^2 k Ak(2jl r jk — 1), 
where 7^ and j k correspond to independent fermionic quasi- 
particles. The energies of these quasi-particles are propor- 
tional to Afc = e' k + A 2 sin 2 k with e k = 1 — A cos fc, and 
the mixing angle appearing in the Bogoliubov transformation 
is 6 k = tan _1 [Asinfc/(e fc + A fe )]. 

The 2D state ^(i)) during the time evolution is the di- 
rect product of N identical ID states \ip(t)). The ID ini- 
tial state \ip(0)) before the quantum quench is the ID ground 
state |o) of the initial Hamiltonian h(0) that is defined by 
Cfe|o) = (Vfc). When this state is expressed in terms of 
the ID ground state of the quench Hamiltonian h(X) that 
is defined by 7fe|w) = (Vfc), it takes the product form 

|V>(0)) - |o) - lL>o( cos ^ " *Bin0fc7£7! fc )|w). Note 
that the connections between the respective 2D and ID ground 
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states are |0) = (g) N |o) and |fi) = 



The state \ip(t)} 



after the time evolution with h(X) is then obtained by the sub- 
stitution 7^, — > -yl_ k e~ 2lAkt in the expression for \ip(0)), and 
the corresponding two-operator expectation values in the po- 
sition representation are given by 



(c4) = m)\ci4\m) = s u . 



{c\cv 



(5) 



{c\c v 



(44: 



TV ^ 



-ik(l-l') s ; n 2 Qk cog 2 Qk gm 2 (2k k t) , 



sin 2 9 k (I 
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k 
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-ik(l-l') 



sin 9 k cos ( 



To calculate the Renyi entropy, we need to evaluate the quasi- 
spin expectation values in Eq. (0). Note that despite the fac- 
torization into independent ID chains in terms of the quasi- 
spins, this calculation gives the Renyi entropy of a strongly 
entangled 2D state |\&(t)) in terms of the physical spins. The 
2D expectation values (^(t)\A x Si A x 2 . . . A x 2ri A x 2r 
are simple products of independent ID expectation val- 
ues (?/.>(i)|A^ Af 2 ...Af Af \ip(t)), where we can write 

A x A x = (4 + c) nfc 1 (! - 2 cJ c i)(4 + c i' ) for each p air ° f 

Af operators. The resulting multi-operator expectation values 
can then be reduced to the two-operator expectation values in 
Eq. §5$ by using Wick's theorem 112211 . 

The long-term time average of the topological Renyi en- 
tropy against the magnetic field A of the quench Hamiltonian 
H{X) is plotted in Fig. [3] The time averages are calculated 
from 1000 independent time instants, which are taken from a 
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FIG. 3: (Color online) Exact long-term time average of 5j as a func- 
tion of/3 = A/(l + A). Three system sizes are plotted: N = 40 and 
D = 3d = 6 (dotted line); N = 60 and D = 3d = 9 (dashed line); 
N = 80 and D = 3d = 12 (solid line). 



sufficiently long interval at sufficiently large times such that 
t 3> max{l, A -1 } for all time instants t. We notice that (Sj) 
is close to the initial value of S-f = 2 for all magnetic fields 
A, and it gets closer to it if we increase the system size TV as 
well as the subsystem dimensions D and d. We therefore ar- 
gue that S2 = 2 in the thermodynamic limit for all times t 
after a quench with an arbitrary magnetic field A. 

Perturbation theory. — In the case of k = 1, the Hamilto- 
nian in Eq. (f2]) requires perturbation theory around the exactly 
solvable point at A = 0. Our aim is to determine the time evo- 
lution of the first order correction to Sj. In general, the most 
straightforward approach is the method of perturbative con- 

~" Q. On the 
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tinuous unitary transformations (PCUT) jl2i Il9l 
other hand, a naive first order perturbation theory is sufficient 
in our case, and we choose to derive our results using this sim- 
pler approach. Note that the same results can also be obtained 
from a PCUT calculation. 

The unperturbed ground state of the initial Hamiltonian 
-ff(O) is |0) and the perturbed ground state of the quench 
Hamiltonian H (A) is |n) as above. In the first order of per- 
turbation theory, the perturbed ground state can be expanded 
as |fi) = |0) + (A/4) £ <ri , r2> |{r a ,r 2 }), where |{n,r a }) is 
the state with two star excitations at the positions ri and r 2 
such that A^ = A Z T = —1. The positions of the stars are 
labeled by the 2D vectors r = (x, y) with < x,y < N. In 
terms of the ground state |n), the initial state |0) becomes 



|*(0)) = |0) = |n>-£ £ |{n,r 2 }) 



(6) 



( r l> r 2> 



Although the non-degenerate corrections to the excited states 
|{ r i; r 2j-) are negligible in the first order calculation, the per- 
turbation introduces a hopping between these degenerate ex- 
cited states. The two star excitations can hop between neigh- 
boring stars with an amplitude — A, and their only interaction 
is a hardcore repulsion: they are not allowed to be at the 
same star. However, since the two excitations are far away 
from each other in the thermodynamic limit, it is reasonable 
to neglect their interaction. Under this assumption, the exact 
eigenstates of the quench Hamiltonian with two star excita- 
tions are |ki,k 2 ) = A r ~ 2 E ri .r 2 e ?;(Ki ri+K2 r2 )|{ri,r 2 }), 
where k = (k x , k y ) with < k x ,k y < N, and each excita- 
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tion has a 2D momentum K = 2n'k/N. Note that the state 
|{ r ii r 2}) = |{i"2,ri}) appears twice in the sum. Since the 
system and the initial state are both invariant under transla- 
tions, we only need to consider the eigenstates with zero total 
momentum ki + k2. These states labeled by |k) = |k, — k) 
have relative energies -E(k) = 4 — 4Ae(k) with respect to 
the ground state where e(k) = cos(K x ) + cos(K y ). In 
terms of the eigenstates, the initial state 1^(0)) takes the form 
|*(0)) = |0) = |0) - (A/2)jXe(k)|k), where the prime 
means that the summation is only over half of the TV 2 possible 
k values. The state after the time evolution with H (A) 

is then obtained by the substitution e(k) —> e(k)e _4iB ^ k ', and 
in terms of the states |{ri , r-2}) it reads 



(7) 



(n- r 2) 



where (ri,r2) means that the summation is over pairs of 
stars without double counting. The coefficient of each state is 
Q(ri,r 2 ,i) = -(A/iV a )X£e(k) cos(K- (n -r 2 )) e - itE W. 
Since the functions Q(ri,r 2 ,i) consist of many incoherent 
oscillations with different frequencies E(k), the average mod- 
ulus square of any such function is given by 



<IO(ri,r 2 ,t)| 2 ) 



cos 2 (K-(n-ra)) = 



A 2 
AN 2 

k 

(8) 

This result has a simple physical interpretation. The two ex- 
citations are always at neighboring stars before the time evo- 
lution, and the sum of the norm squares in the states with two 
excitations is 2iV 2 A 2 /16 = N 2 X 2 /8 because there are 2N 2 
ways of choosing two neighboring stars. During the time evo- 
lution, the excitations hop between stars, and the norm square 
is distributed uniformly between all possible states with two 
star excitations. Since there are N 2 (N 2 - l)/2« N 4 /2 ways 
of choosing any two stars, the uniform norm square corre- 
sponds to the result X 2 /AN 2 in Eq. ©. 

When calculating the Renyi entropy for the state \^(t)} in 
Eq. (0, there are two corrections to the value 5f A = L—n for 
the unperturbed ground state |0): a static correction from the 
perturbed ground state | f2) and a dynamic correction from the 
oscillating terms oc Q(ri, r 2 , 4). The static correction is lin- 
early proportional to the combined length L of the boundary in 
each case (m), therefore its topological contribution vanishes 
1 1211 . On the other hand, the terms oc Q(ri, r 2 , t) give an ex- 
pectation value (*(t)|i^ij 2 |*(t)) = 2Re[Q(ri,r 2 ,f)] for 
each pair of stars ri and r 2 on the same closed loop of the 
boundary. Since (Re[Q(n, r 2 , t)] 2 ) = (|Q(n, r 2 , t)\ 2 ) /2 
and there are £{£ — l)/2 ways of choosing two stars from a 
closed loop of length I, the time average of the dynamic cor- 
rection is (ASf A ) = -X 2 £(£ - 1) /{AN 2 In 2) for each closed 
loop of the boundary in each case (m). Since there is a term 
oc £ 2 , the topological contribution does not vanish. In particu- 
lar, when the subsystem dimensions satisfy D = 3d, the time 
average of the topological Renyi entropy becomes 

T SGA 2 ^ 2 
^ 2 >- 2 7V2 ln2 - (y) 




FIG. 4: (Color online) Long-term time average of S2 from the ex- 
act treatment (solid line) and from the first order perturbation the- 
ory (dashed line). Two system sizes are plotted: TV = 40 and 
D = 3d = 6 (a); N = 80 and D = 3d = 12 (b). 



If we take both the system size N and the subsystem dimen- 
sion d to infinity such that d ~ \fN, the perturbative cor- 
rection vanishes as 1/N. We can thus choose macroscopic 
subsystems in the thermodynamic limit in such a way that the 
topological Renyi entropy is 5j = 2 for all times t. 

To verify that the first order perturbation theory is indeed 
a good approximation at A <C 1 and it is not invalidated by 
higher order corrections, we also establish an analogous per- 
turbation theory in the exactly solvable k — case. When 
D = 3d, the time average of the topological Renyi entropy 
after the first order correction takes the form 



(.5': 



A 2 (d 2 + 12d- 1) 



2iVln2 



(10) 



The exact results and those obtained from the perturbation the- 
ory are plotted together in Fig. |U We notice that the respective 
curves are in good agreement when A < 1, and therefore we 
argue that the first order perturbation theory is a good approx- 
imation in the k = 1 case as well. 

Conclusions. — In this Letter, we studied the resilience of 
topological order in the TCM ground state after a quantum 
quench with an external magnetic field. The time evolution 
of topological order was detected via the topological Renyi 
entropy of order 2. We considered two different quenches: the 
integrable one was solved exactly, while the non-integrable 
one was treated with perturbation theory. In both cases, we 
found that topological order is resilient in the thermodynamic 
limit. Our results are in agreement with those in Ref. 11811 . but 
they apply to significantly larger system sizes. 

It is interesting to discuss the generic conditions under 
which topological order is resilient. First of all, the results 
for the non-integrable quench show that integrability does not 
play a role. On the other hand, both quenches preserve the Z2 
gauge structure of the TCM ground state. This gauge struc- 
ture allows us to use the modified definition for S%, which in 
turn makes the subsequent calculations possible. We therefore 
believe that topological order characterized by Sj is a robust 
non-equilibrium property of the system whenever gauge in- 
variance is preserved. In perspective, it would be important to 
study the behavior of topological order after a more generic 
non-gauge-preserving quench. This further step could help us 
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better understand the universal non-equilibrium and thermal 
properties of topologically ordered systems. 

Acknowledgments. — We thank X.-G. Wen for illuminating 
discussions. This work was supported in part by the National 
Basic Research Program of China Grant (201 1CBA00300, 
2011CBA00301) and the National Natural Science Founda- 
tion of China Grant (61073174, 61033001, 61061130540). 
Research at the Perimeter Institute for Theoretical Physics 
is supported in part by the Government of Canada through 
NSERC and by the Province of Ontario through MRI. 



[1] X.-G. Wen, Quantum Field Theory of Many-body Systems (Ox- 
ford University Press, 2004). 
[2] A. Y. Kitaev, Ann. Phys. (N. Y.) 303, 2 (2003). 
[3] M. H. Freedman, A. Kitaev, and Z. Wang, Commun. Math. 

Phys. 227, 587 (2002); C. Nayak, S. H. Simon, A. Stern, M. 

Freedman, and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008). 
[4] A. Hamma, R. Ionicioiu, and P. Zanardi, Phys. Lett. A 337, 22 

(2005); A. Hamma, R. Ionicioiu, and P. Zanardi, Phys. Rev. A 

71, 022315 (2005). 
[5] A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 (2006); 

M. Levin and X.-G. Wen, ibid. 96, 1 10405 (2006). 
[6] S. T. Flammia, A. Hamma, T. L. Hughes, and X.-G. Wen, Phys. 

Rev. Lett. 103, 261601 (2009 ). 
[7] I. H. Kim JarXiv: 1210.23601 

[8] A. Hamma, W. Zhang, S. Haas, and D. A. Lidar, Phys. Rev. B 

77, 155111 (2008). 
[9] S. V. Isakov, M. B. Hastings, and R. G. Melko, Nature Phys. 7, 



772 (2011). 

[10] L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. 

Phys. 80, 517 (2008). 

[11] P. Zanardi and L. C. Venuti. larXiv: 1205.25071 

[12] G. B. Halasz and A. Hamma, larXiv: 1 206.4223 1 

[13] P. Horodecki and A. Ekert, Phys. Rev. Lett. 89, 127902 (2002); 

S. J. van Enk and C. W. J. Beenakker, ibid. 108, 1 10503 (2012). 
[14] D. A. Abanin and E. Dernier, Phys. Rev. Lett. 109, 020504 

(2012). 

[15] A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, 

Rev. Mod. Phys. 83, 863 (2011). 
[16] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and I. 

Bloch, Nature 415, 39 (2002); M. Greiner, O. Mandel, T. W. 

Hansch, and I. Bloch, Nature 419, 51 (2002). 
[17] P. Calabrese and J. Cardy, Phys. Rev. Lett. 96, 136801 (2006); 

M. A. Cazalilla, ibid. 97, 156403 (2006); M. Rigol, V. Dunjko, 

V. Yurovsky, and M. Olshanii, ibid. 98, 050405 (2007). 
[18] D. I. Tsomokos, A. Hamma, W. Zhang, S. Haas, and R. Fazio, 

Phys. Rev. A 80, 060302(R) (2009). 
[19] S. Dusuel, M. Kamfor, K. P. Schmidt, R. Thomale, and J. Vidal, 

Phys. Rev. B 81, 064412 (2010). 
[20] I. Yu, S.-P. Kou, and X.-G. Wen, Europhys. Lett. 84, 17004 

(2008). 

[21] E. Barouch, B. M. McCoy, and M. Dresden, Phys. Rev. A 2, 
1075 (1970). 

[22] E. Barouch and B. M. McCoy, Phys. Rev. A 3, 786 (1971). 

[23] F. Wegner, Ann. Phys. (Leipzig) 3, 77 (1994); S. D. Glazek and 
K. G. Wilson, Phys. Rev. D 48, 5863 (1993); S. D. Glazek and 
K. G. Wilson, ibid. 49, 4214 (1994); S. Dusuel and I. Vidal, 
Phys. Rev. B 71, 224420 (2005); S. Dusuel, M. Kamfor, R. 
Onis, K. P. Schmidt, and J. Vidal, Phys. Rev. Lett. 106, 107203 
(2011). 



